Exact dynamics in the inhomogeneous central-spin model 
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We study the dynamics of a single spin- 1/2 coupled to a bath of spins- 1/2 by inhomogeneous 
Heisenberg couplings including a central magnetic field. This central-spin model describes decoher- 
ence in quantum bit systems. An exact formula for the dynamics of the central spin is presented, 
based on the Bethe ansatz. For initially completely polarized bath spins and small magnetic field 
we find persistent oscillations of the central spin about a nonzero mean value. For a large number of 
bath spins Nb, the oscillation frequency is proportional to Nb, whereas the amplitude behaves like 
1/Nb, to leading order. No asymptotic decay of the oscillations due to the non-uniform couplings is 
observed, in contrast to some recent studies. 

PACS numbers: 73.21.La, 03.65.Vf, 02.30.1k 
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I. INTRODUCTION 

The central-spin model, or Gaudin model 1,2 , is defined 
through the Hamiltonian 



;V(, 



(1) 



One central spin So interacts with Nb bath spins via in- 
homogeneous Heisenberg couplings Aj. Additionally, a 
magnetic field h is included that couples to the central 
spin only. Here, we restrict ourselves to spin-1/2 objects. 
This model is being used widely to model the hyper- 
fine interaction between an electron spin in a quantum 
dot and the surrounding nuclear spins 3-11 . The localiza- 
tion of the electron wave function within the quantum 
dot leads to inhomogeneous coupling constants, depend- 
ing on the distance of the nuclei from the center of the 
dot ' 8 . Since electrons in quantum dots are candidates 
for spin qubits 12 , special interest is on decoherence phe- 
nomena of the central spin. In ,s it was argued that the 
non-uniform hyperfine couplings cause the central spin 
to decay. Two special cases were considered in those 
works: a completely polarized and a completely unpo- 
larized bath. In the case of a strong /i-field, results for 
arbitrary polarization have been obtained in 13 . In all 
cases, the initial oscillations of the central spin were re- 
ported to decay, implying that the central spin comes to 
rest at some equilibrium value due to the inhomogeneous 
couplings to the environmental spins. In this paper, we 
show that at least for an initially completely polarized 
bath, the exact solution from the Bethe ansatz (BA) be- 
haves differently. 

The model Eq. (1) has also long been of interest from 
a fundamental point of view. As mentioned above, its 
eigenvalues and eigenstates can be constructed exactly 
in the framework of the BA 1 ' 2 ' 14 ' 15 . In this formalism, 
the solution of the many-particle eigenvalue problem is 



achieved by the solution of coupled algebraic equations 
for the quantum numbers (BA numbers). Both the eigen- 
states and the eigenvalues are constructed in terms of 
these numbers. However, no BA studies of the dynam- 
ics have been performed yet, which is mainly due to the 
intricate distribution of the BA numbers in the complex 
plane. 

Besides the full quantum-mechanical solution, a com- 
plementary approach consists in treating the model (1) 
on a mean-field level, where operators arc replaced 
by their expectation values. It was argued 11,16 that 
this approach is exact in the thermodynamic limit. 
But, as pointed out very recently 17 , agreement between 
the quantum and the mean-field solutions depends on 
the choice of the initial conditions in the quantum- 
mechanical solution. In the question was addressed to 
what extent initial entanglement between the bath spins 
is essential for the mean-field behaviour. If all couplings 
are equal, Aj =: A\/j in (1), then in the extreme case 
of an initially unentangled unpolarized bath the central 
spin was observed to decay from its initial value ±1/2 to 
a value ~ 1/AT& after a decoherence time ~ l/^/TV^. 
In other words, in the thermodynamic limit, decoher- 
ence occurs infinitely fast and after that, the dynamics is 
frozen. The aim of the present paper is to shed more light 
on the time evolution starting from an unentangled bath 
state for the physically relevant case of inhomogeneous 
couplings, using the exact BA solution. 

During the past few years, a whole variety of meth- 
ods has been applied to study decoherence phenomena 
in the central-spin model. However, the potential of the 
BA method is still unexplored. In this regard, this work 
constitutes a first step to demonstrate the practicability 
of the BA approach by first presenting a general formula 
for the dynamics of the central spin and then focusing on 
one special case. The study of this special case should 
be considered not as the end, but rather as the starting 
point of the investigation of the exact solution. 
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In the following, we derive an exact formula for the 
time-dependent expectation value (S§)(t) of the central 
spin by employing the BA solution. That formula is our 
central result. It is valid for an initial state which is a 
mutual eigenstate of all 5| operators and hence a prod- 
uct state containing no entanglement at all. The degree 
of spin polarization of the initial bath state is an essential 
control parameter. For M b flipped spins (as compared to 
the ferromagnetic "all up" state) in the bath and with 
the central spin flipped as well, the result for (Sq)^) is 
expressed in terms of a sum of harmonic oscillations with 
frequencies determined by the sum of M b + 1 BA num- 
bers, and with amplitudes given through the norms of the 
eigenfunctions, see Eqs. (6,7) below. We evaluate (Sq)^) 
analytically for M b = (fully polarized bath), with cou- 
plings distributed in an interval ]0,„4] (A =const). Un- 
der these conditions, (S^HO oscillates with a frequency 
~ \ Y^,j=i Aj and an amplitude of order 0(1/N b ) about 
a mean value (Sg){t) = -± + 0(1/N b ). The ampli- 
tude initially decreases but stays constant after a time 
t a ~ Air/ A. This differs from the results of 7,8 , where the 
oscillations of (Sq)^) were found to decay completely 
Our analytical evaluation of (S§)(t) is confirmed by a 
numerical evaluation for zero field and N b = 30 bath 
spins. 

An analytical evaluation is also possible for a large 
magnetic field (jh\ 3> N b ), with similar results, namely, 
persistent oscillations of amplitude 0(N b /h 2 ). From the 
magnetic field dependence of the results we conclude that 
there exists a resonant value for h, where (S*o)(i) oscil- 
lates at the maximum amplitude possible. This result is 
also confirmed by numerical evaluation for N b = 30. 

The remainder of this paper is organized as follows. In 
the next section, we sketch the BA setting and give the 
general formula for (Sq)^). The third section contains 
both analytical and numerical results for a fully polarized 
bath. 



II. THE BETHE ANSATZ SOLUTION 

The Bethe ansatz solution of the model (1) was found 
by Gaudin 1,2,15 . Before summarizing it, let us introduce 
the abbreviations N b A := Y^fli Aj, N := N b + 1, M := 
M b + 1, C M := N\/(M\(N-M)\). The eigenvalues k„ of 
the Hamiltonian (1) for a fixed number M of down-spins 
(i.e. for magnetization N/2 — M) were given in 2 ' 14 in 
terms of the BA numbers w^y? 1 

1 M N h 

fc=l 

provided that the tu kL , fulfill the equations 
Nh , m 



An 



J Aj — ujk,u 



2h 



= 



for k = 1, . . . , M. Gaudin 15 showed that there are C M 
sets of solutions . . . , ujm.v} to these equations, one 

for each eigenvalue A„. The corresponding energy eigen- 
states with a fixed number M of flipped spins read 



M N b 
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\ ^k,v - A jk 



S~\0), (3) 



where Aq = oo by definition. Before writing down 
the normalization factor tim v , let us comment on the 
last formula. By |0) we mean the fully polarized state 
| ft;T, . . ■ ,t), where the symbols ft, Jj. for the central spin 
and |, J. for the bath spins are used. In the second 
step of (3), the product of the previous equation is ex- 
panded into a sum over different ordered configurations 
J = {ji, . . . , Ja/}, where M spins are flipped at the sites 
ji,...,M- The first sum in the final expression of (3) com- 
prises the M\ permutations of the set J. The normal- 
ization factor tim v was conjectured by Gaudin 2,15 and 
proved by Sklyanin 18 for h = 0: 
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For M = 1 this is obviously true also for finite h. We 
have checked the validity for N b = 3, M = 2 and finite 
h as well and thus conjecture that this formula holds for 
general N b , M, h. 

In order to arrive at an expression for (Sq)^), we de- 
compose the initial state into Bethe ansatz eigenstates 
(3). We focus here on the case where the initial state 
is a product state. It is denoted by \L), where L = 
{£i, . . . ,1m} is the set of lattice sites with initially flipped 
spins. Let us rewrite Eq. (3) by introducing the C M x C M 
matrix D such that 



Ml M 



\M V ) =Y^D v ,j\J) , D v ,j = — e n 
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where J again denotes the set of locations of flipped spins. 
From the hcrmiticity of the Hamiltonian and from the 
normalization of the eigenstates it follows that D is uni- 
tary: [-D -1 ]^ j = Dj v , where * denotes complex conju- 
gation. Thus the time evolution of the initial state reads 



\L)(t) = 5>; L A,./e- iA "V>- 



(4) 



We assume that the central spin is initially flipped, i.e. 
Mb = M — 1. Then l\ = and the initial configuration 
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is completely determined by the set Lb = {£2, ■ ■ ■ ,@m} 
which contains bath sites only. Ordered sets Jf, can be 
defined analogously. From Eq. (4), one can infer the re- 
duced density matrix for the central spin 



Po (t) = $> j6 (0| 2 m<4 



with 



7 , A/,0UZ, t ,- t> V,0UJi, 



D„ 



(5) 



(6) 



The ordered set U contains the central spin site 
and Mf, bath sites; U Jf, is defined similarly. We also 
inserted the eigenvalue (2) and dropped an overall phase 
factor in the last equation. From Eq. (5), we conclude 
that 



/ 



(7) 



\ 



Let us pause briefly to comment on the structure of (7) 
which is the central result of this paper. Each aj b (t) 
is a superposition of exponentials e _1 * A " , with each 
Aj, the sum of M = Mb + 1 Bethe ansatz numbers uJk,v, 
apart from trivial constants and factors, see (2). The ex- 
pectation value (7) thus contains oscillations with combi- 
nations of all these frequencies. It depends on the initial 
state through the dependence of |aj 6 (t)| in Eq. (6) on 
Lb, the set of initially flipped bath spins. Let us also 
comment on the Poincare recurrence time rp after which 
the system returns to the initial state. Generally, BA 
numbers uik, v are complex (either real or complex con- 
jugate pairs), such that the eigenvalues A„, Eq. (2), are 
irrational numbers. This means that the system never 
reaches its initial configuration again, i.e. rp — > 00. In 
special cases, however, the recurrence time can be finite, 
for example for homogeneous couplings, Aj = A 17 , where 
an exact solution without the BA is possible. Thus in 
that particular case, the BA numbers must be rational, 
which we verified for Mb = 0, N arbitrary and for Nb = 3, 
M b = 1. 



III. SPECIAL CASES 

In the homogeneous case 22 (Aj = AV j) the Hamilto- 
nian can be diagonalized directly, i.e. without employing 
the Bethe ansatz 17 . We have checked that for Nb arbi- 
trary, M b = and for N b = 3, M b = 1, Eq. (7) yields the 
same results as those obtained in 17 . 



A. Fully polarized bath: Analytical results 

Let us now concentrate on the inhomogeneous case 
with arbitrary Nb and a fully polarized bath (Mb = 0). 
The initial state then is | JJ-; f ... f). The sum in Eq. (6) 
simplifies to 



N 



«(*) = E 



e W, 



(8) 



where the frequencies oj v are solutions of the single BA 
equation 23 



A, 



y 



2h 



and (<So)(t) is given by 
(W) = ; 



(l-2|a(i)| s 



(9) 



(10) 



Before solving these equations numerically, we first ex- 
tract the qualitative behaviour of (S§)(t) in the limit 
of large Nb. Therefore, the N roots uj v of (9) have to 
be found. In all that follows, we restrict ourselves to 
the physically relevant case Aj > OVj and assume the 
couplings to be distributed such that for large particle 
numbers, N b A - N b and £ M - N b . 



Weak field 



Let us first consider a weak magnetic field, \h\/N <C 



infjAj. Then 



uji = -2h/N + 0(h 2 /N 2 ) 



111) 



is one solution, tending to as h — > 0. Another solution 
is found at 



uj 2 = N b A-2h + 0(l) 



(12) 



By sketching a graphical solution of Eq. (9) one sees that 
each of the remaining TV — 2 solutions is located between 
two consecutive Aj. 

For the sake of simplicity, we write down a(t) first for 
h = 0, and calculate the amplitudes in (8) only to order 
0(l/N b ). 



a(t) 



A 



N 

E 

v=3 



(N b A) 2 

N,, 



1 e ±N b At 



i + E 



4 



(Aj 



1 



(13) 



+0 (1//V 2 ; 



To obtain a crude approximation of the last sum, we 
recall that each oj„ is located between two consecutive An . 
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Let < Aj < A, where A is independent of N, and let 
the Aj be distributed smoothly between and A. Then 
we estimate the difference between two consecutive Aj as 
~ A/Nb, such that those Aj which are closest to u>„ in 
the last bracket in (13) dominate and yield a contribution 
~ A 2 /N 2 to each term of the sum over the u> v . We 
assume here that the number of these terms does not 
scale with Nb such that the whole sum behaves as 



N b 
4=1 ^ 3 



42 



N 



b ■ 



(14) 



Obviously, we cannot perform the thermodynamic limit 
here directly since then the sum (14) would diverge. Con- 
sequently the third term in (13) would vanish, as would 
the second term, resulting in trivial dynamics. (Note that 
|a(t)| 2 would equal unity in that case.) Thus a(t) above 
has to be evaluated for large, but still finite Nb, avoiding 
singularities that occur for N b — > 00. 

Therefore the number of particles is kept finite, and, 

r jy ^2 

following Eq. (14), we replace 1 + ^ h 
IvjN 2 , with certain constants 7„ 



Then 



(Aj- 



N 

E 

i/=3 



N, 



1 N 



(15) 



We will now show that this sum describes an amplitude- 
modulated oscillation where the envelope has its first zero 
around 



t a := 4tt I A 



(16) 



For a large number of particles, we approximate the last 
sum in Eq. (15) by an integral. To leading order this 
yields 

N f A 

Y,J^ ,A ~N b / /(w)e^p(o;)dw + O(l), 

where f(tu) is the continuum limit of the 7^ in Eq. (15) 
and p(lu) is the root density. Proceeding further, we can 
write 

JV 



u=3 



r A/2 



e iHA)t I ■ + A/2)e^ ult p(Lu + A/2)du; 

J -A/2 



(17) 



where the constant b(A) in the phase factor is chosen 
such that Fa is real (if f(to) and p(uj) are constants, then 
b(A) = -4/4). The integral is the Fourier transform of a 



band-limited function of bandwidth A. The frequency 
cutoff at ±A/2 leads to a decaying oscillation with a 
period 8tt/A, which reaches its first zero around t a . By 
inserting the result (17) into Eq. (13), we evaluate |a(£)| 2 
including orders 1, l/Nt,: 



\a(t)\' 



(N b A)* ^ 



E4 



(18) 



2 



cos 



N h A 



Fj^(t) cos 



N h A 



- b(A) t 



Thus the qualitative behaviour of |a(i)| 2 is as follows. 
At t = 0, the initial condition is |a(0)| 2 = 1. After that, 
|a(t)| 2 displays fast oscillations around the mean value 



E^: 



{N b AYf^ 
with the leading frequency 



NhA 



(19) 



(20) 



These oscillations are given by the two cosine terms in 
Eq. (18). Both terms have amplitudes of order 0(l/N b ). 
But whereas the amplitude of the first cos-term is con- 
stant, the amplitude of the second one depends on time. 
Especially, if ^vi(O) > 1, the latter term will dominate 
initially, decaying ~ t 2 for short times t <^t a . As stated 
above, J-A(t) decays and it is certainly dominated by the 
other cosine term at times t ~ t a . As the exact numerical 
evaluation below (see Fig. 2) shows for both a uniform 
and a nonuniform distribution of the couplings, the con- 
stant amplitude term also dominates for all times t > t a . 

In any case, |a(i)| 2 does not show any long-term decay, 
and consequently neither docs (S^)(t). We rather expect 



Ht)f 



{N b AY ^ 



NhA 



h t 



(21) 



for long times t^$>t a , with an undetermined amplitude c\ 
of order 0(1). Although the calculations leading to the 
results (18,21) certainly lack mathematical rigor, they are 
confirmed in the next section by a numerical evaluation 
of the exact Eqs. (8, 9). 

Note that the last sum in Eq. (13) is completely due 
to the inhomogeneities of the couplings: If all couplings 
are equal, the Bethe Ansatz equation (9) is a quadratic 
equation in u> with the two solutions u>i t 2 (the exact so- 
lutions in that case are given below). Summarizing, the 
non-uniformity of the couplings causes an initial decay of 
the amplitude until a time ~ t a . For longer times, the 
amplitude does not change any more and (So)(i) keeps 
oscillating with an amplitude of order 0(1/N b ) about a 
mean value — \ + 0(1/N b ). This docs not agree with 7,8 , 
where the amplitude of oscillations is predicted to vanish 
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in the long-time limit. The leading frequency of the oscil- 
lation is uji (20), in agreement with 7,8 . A small positive 
magnetic field causes a slow-down of the oscillation, as 
can be seen from Eq. (12). 



2. Strong field 



with 



A 



A 2 & 



(A+A 2 /(2h))t~ 



p{A)dA. 



We now discuss the case of a strong field \h\ ^> N b . 
Here the limiting values of the two solutions oj\^ consid- 
ered above in the weak field case, Eqs. (11,12), are given 
by 



wi = -2h + N b A + 0{N b /h) 
uo 2 = A l (l + A l /{2h)) + 0(l/h 2 ). 

The remaining BA roots are found at 

ccv = A„_i(l + A u - 1 /{2h)) + 0{l/h 2 ), u = 3,. 



(22) 
(23) 



A(.24) 



Thus the magnetic field shifts the roots along the real 
axis. (See Fig. 1 for an example.) 

Eqs. (22,23,24) are inserted into Eq. (8), keeping only 
the leading non-trivial contribution in the amplitudes: 



a(t) 



1 



- w ,Y. A ^^ {h - NbA/2)t 



^ Ah 2 



e ^(A J +A 2 J /(2h))t^ 



Arguing as in the zero-field case after Eq. (13), we find 

N b 

2h 2 



\a(t)f 



1 



1 N b 



2Nh r 
' h 2 



F%\t)\ cos [(h-N b A/2 + b h {A)) t] ,(25) 



Here p(A) is the density of the coupling constants A/-, and 
bh{A) in Eq. (25) is chosen such that Fj^ (t) is real. Note 
that in this case, the result depends only on the (known) 
couplings Ak . It is expected to be valid for times as long 
as the lo v can be approximated by Eqs. (22,23,24), i.e., 

t <C h/Nb- The enveloping function Fj^'(t) again de- 
scribes a decreasing oscillation. For times t 3> h/N b , the 
replacement of the ui v by the A v -\ is no longer valid, and 
|a(i)| 2 oscillates around a mean value 1 — ^2 Ylj=i 
with an amplitude ~ N b /h 2 and frequency ~ \h—N b A/2\. 
Note that increasing a strong field suppresses the ampli- 
tude further and raises the frequency. Combining this 
with the findings for the small field case, we conclude 
that for h > 0, there should be a resonance region around 
a field h r where the amplitude is maximal and the fre- 
quency minimal. 

3. Resonance field 



To study that region, it is instructive to consider ho- 
mogeneous couplings first, Aj = AVj. In that case the 
BA equation (9) becomes a quadratic equation with the 
two solutions u>\ t 2 given below. The remaining solutions 
^v=3,,....n all tend to u> v = A when Aj — > A\/j, as can 
be shown using techniques from 15 . The contributions to 
a(t), Eq. (8), corresponding to those roots all vanish and 
we end up with 



Ht)\ 2 

Wl,2 



-4 



LOl — LU2 



^ (AN -2hT V 8Ah + ( AN ~ 2h ) 2 



1 



2 ((AN - 2h) 2 + 8Ah) 



(AN - 2(A + h)) 2 + AA 2 (N - 1) cos 



cJi - 0J2) t/2] 



-y/(AN - 2h) 2 +8Aht 



For h = A(N - 2)/2 =: h r , this results in S§(t) = 
— I cos [Ay 7 ^ t] and 

^u=4(-1tv^)- (26) 

Thus if h equals the C(7Vf,)-magnctization of the bath, 
the central spin resonates with maximal amplitude and 



the two roots uj\,2 have equal absolute values in the limit 
of a large particle number. The width of the resonance is 
estimated by considering the time averaged value (Sq) as 
a function of h. From this one deduces a width ~ y/Nb- 
This behaviour was also found in 7 ' 8 for the inhomogc- 
neous case, which we turn to now. 
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By analogy to the homogeneous case we expect reso- 
nance to occur at 



h r = N b A/2. 



(27) 



^From the discussion of the large- and small-/i limits 
above we expect that at least one \u) u \ will be much larger 
than all Aj. Expanding the BA equation (9) to second 
order in lo~ x we obtain a quadratic equation for the ex- 
tremal u>„: 



E =1 A j 



L,j=i A 



2h 

<jj„ 



(28) 



For the expected resonance field value (27) this leads to 



"1,2 



0(1). 



(29) 



Hence the two extremal BA roots W1.2, for which limit- 
ing values for small and large fields have been given in 
Eqs. (11,12) and (22,23), respectively, have equal abso- 
lute values, and the resonant behaviour is determined (to 
leading order) by 



|a(t)| 5 



1 + cos (cJi,2i)] 



(30) 



The width of the resonance is still ~ \ Nb, in agreement 
with 7 ' 8 . 

We illustrate the evolution of the BA roots with h by 
solving Eq. (9) numerically for Nb = 10 and Aj = (11 — 
j ) / 1 , j = 1, .. . , 10 and h between and 5. The results 
are depicted in Fig. 1. As expected, the "inner roots" 
w 3,...,w are always located between two successive Aj, 
whereas the location of the "outer roots" u>i t 2 changes 
significantly with h. 



B. Fully polarized bath: Numerical results 

In the following, we show exact results for (Sq) (t) with 
Ni, = 30 bath spins, obtained by numerically solving the 
BA equation (9), and compare those results with the ana- 
lytic approximations just derived. Note that the features 
of those approximations are universal in the sense that 
they depend on the distribution of the Aj only through 
the mean value A. We consider two different choices of 
distributing the Aj between zero and A = 1; a nonuni- 
form distribution and a uniform one: 

Aj = cxp[-(j/a) 2 ] witha = iV&/2 (31) 
Aj = (N b + l-j)/N b (32) 

for j = l,...,N b . 

Figure 2 shows the time evolution of (Sq) for zero field, 
h = 0. The initial decay ~ t 2 of the envelope, the time 
scale t a where that decay stops, the frequency and the 
mean value around which (<5o)(X) oscillates agree well 





1 


' 1 


4 






d 

-2 












-4 








1,1 


iii, 



FIG. 1: The locations of the N Bethe ansatz numbers u v for 
N = 11 and a uniform distribution of the couplings between 
and 1. The (blue) straight diagonal lines to the left and right 
are the small- and large-field approximations Eqs. (12,22). 
The (red) vertical line denotes the resonance field h r , Eq. (27), 
and the two crosses are the approximate values for wi,2 from 
Eq. (29). 



with the exact numerics. We recall that from combin- 
ing Eqs. (10,19), the mean value follows 



N b 



{NbAf^ 



(33) 



whereas the leading (i.e. the fast) oscillation is given by 
LOi, Eq. (20). In order to compare that (approximate) 
value with the exact numerics we counted the oscilla- 
tions within the time interval shown. The numerical val- 
ues thus determined are u>i :n = 6.7, (w; >n = 8.1) for the 
nonuniform (uniform) distributions, in good agreement 
with the analytical predictions loi = = 6.4 (7.8). 

The resonance case h = h r is shown in Fig. 3. The os- 
cillation is between ±1/2, as expected from Eqs. (10,30). 
The numerical values for the frequencies uj^ n = 3.0 
(u>l t n = 3.3) for the nonuniform (uniform) distribution 
compare well with the approximate values derived from 

Eqs. (29,30), ySS^f = 3.0 (3.2). 

For the case of a large field, we checked that the ini- 
tial decay is ~ t 2 , in agreement with Eq. (25). Fur- 
thermore, the approximate mean value (Sq) = —1/2(1 — 
pr J2f=i A ]) coincides with numerical results as well. 



IV. CONCLUSION AND OUTLOOK 

The approximate analytical and exact numerical eval- 
uations of the general formula for the time evolution of 
the central spin polarization {Sq) for an initially fully 
polarized spin bath show that in this case, an inhomo- 
geneous broadening of the Heisenberg couplings leads to 
decohcrence only initially for short times. As exemplified 
by Fig. 2, this decohcrence process is far from complete 
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FIG. 2: (So)(t) for h = 0, N b = 30 bath spins and a non- 
uniform (upper panel) and uniform distribution (lower panel) 
according to Eqs. (31), (32) respectively. Red dashed lines 
denote the mean value (33), and blue dotted lines the time 
scale t a , defined in Eq. (16). The initial decay of the envelope 
is ~ t 2 , as expected from Eq. (18). 



-0.25 
-0.5 



and does not suppress the oscillations of the central spin 
completely at long times. 

Our previous work 17 on the central spin model with ho- 
mogeneous couplings shows that an initially uncntanglcd 
spin bath supports complete decohercncc, if the magne- 
tization of the bath is zero or small. In that case (S§)(t) 
decays to zero within a decoherence time ~ 1/ \/N b . 
However, at a later time rp = 0(1) (the Poincare re- 
currence time) (S§)(t) shows a complete revival due to 
the commensurate energy spectrum of the homogeneous 
model. In the inhomogeneous model we expect a diver- 
gent Poincare recurrence time, Tp — > oo. 

In contrast to that an entangled bath state with zero 
or small magnetization may lead to persistent oscillations 
of (So)(t) at maximum amplitude in the homogeneous 
model. We thus see that in the homogeneous model ini- 
tial bath states with zero magnetization but different de- 
grees of entanglement may lead to very different long- 
time behaviors. This observation is consistent with 19,20 , 
where it was argued that entanglement in the bath pro- 
tects the central spin from decohering. If this scenario is 
generally valid remains to be studied further. 

Another very interesting question is to what extent 
non-uniformity of the couplings affects the decoherence 
of the central spin for an unmagnetized bath, as opposed 
to the completely magnetized bath studied in the present 
paper. Although this question has been addressed with 
perturbative methods before 7,8 , the analysis of the gen- 
eral solution derived in the present paper is expected to 
yield deeper and more quantitative insights. 



FIG. 3: {S§)(t) for h = h r , N b = 30 bath spins (cf. Eq. (27)) 
and a nonuniform (upper panel) and uniform distribution 
(lower panel) according to Eqs. (31), (32) respectively. In this 
case, h r = 6.37 for the nonuniform and h r = 7.75 for the 
uniform distribution. 
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It is worthwhile noting that these equations can also be 
obtained by the ansatz made in 7 ' 8 . In that approach, the 
BA numbers are the poles of the Laplace transform of a(t). 
In contrast to 7,8 , where the many-particle limit is taken in 
the Laplace domain, here this limit is taken in the time- 
domain. 



